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Abstract 

We consider option pricing in a regime-switching diffusion market. As 
the market is incomplete, there is no unique price for a derivative. We 
apply the good-deal pricing bounds idea to obtain ranges for the price of a 
derivative. As an illustration, we calculate the good-deal pricing bounds 
for a European call option and we also examine the stability of these 
bounds when we change the generator of the Markov chain which drives 
the regime-switching. We find that the pricing bounds depend strongly 
on the choice of the generator. 



1 Introduction 

Regime-switching market models are a way of capturing discrete shifts in mar- 
ket behavior. These shifts could be due to a variety of reasons, such as changes 
in market regulations, government pohcies or investor sentiment. First intro- 



duced by Hamilton| [1989j , regime-switching models have been shown in various 



empirical studies to be better at capturing market behavior than their non- 



regime-switching counterparts (for example, see Ang and Bakaert 2002 , Gray 
[T996j and Klaassen 2002j ). 

An example of regime-switching market is one in which there are only two 
regimes: a bear market regime and a bull market regime. Suppose the market 
starts in a bull market regime, in which prices are generally rising. It stays 
in this regime for a random length of time before switching to a bear market 
regime, in which prices are generally falling. It then stays in the bear market for 
another random length of time before switching back to the bull market. This 
cycle continues ad infinitum. 

Due to the regime-switching, the market is incomplete and hence there is no 
unique risk-neutral martingale measure to use for pricing derivatives. In fact, 
there are infinitely many possible risk-neutral martingale measures. This means 
that not only is there no unique price for derivatives, but the range of prices 
which can be obtained from all the possible risk-neutral martingale measures 
are too wide to be useful in practice. 

As prices of derivatives are not unique in incomplete markets, various sug- 
gestions have been made either on how to choose a single price or on how to 



obtain a more restricted, and therefore potentially more useful, range of prices. 
We focus in this paper on the latter because it is the market which ultimately 
decides which risk-neutral martingale measure is used for pricing a derivative 
and so we should take into account our uncertainty about what the market 
price will be. Therefore, we believe it is better to find a range of prices that 
the market-determined price might reasonably be expected to lie in, rather than 
determining a single price. 

The idea that we build upon is that of the good- deal bound. This idea is 
due to Cochrane and Saa Requejo 2000 and is based on the Sharpe Ratio, 



which is the excess return on an investment per unit of risk. Their idea is to 
bound the Sharpe Ratios of all possible assets in the market and thus exclude 
Sharpe Ratios which are considered to be too large. The bound is called a 
good-deal bound. The method of applying the bound gives a set of risk-neutral 
martingale measures which can be used to price options. This results in an upper 
and lower good-deal bound on the prices of an option. The idea was streamlined 
and extended to models with jumps in Bjork and Slinko 2006 , and it is their 
approach that we follow in this paper. 

The uncertainty measured by the size of the good-deal bounds reflects the 
uncertainty within the market model concerning the price of the derivative. It 
does not measure uncertainty concerning the choice of the model, by which we 
mean both the model structure - in this case a regime-switching model - and 
the model parameters. Indeed, as we see concretely in a numerical example, 
changing the model parameters changes a derivative's pricing bounds for a fixed 
good-deal bound. Thus the model choice is still a very important factor in 
determining the derivative's good-deal pricing bounds. In summary, the good- 
deal pricing bounds tell us nothing about model uncertainty itself, only about 
the uncertainty in the choice of a risk-neutral martingale measure within a 
particular model. 

outline various ways that the good-deal 



Cochrane and Saa Requejo 2000 



pricing bounds can be used, such as a trader using the bounds as buy and sell 
points and a bank using them as bid and ask prices for non-traded assets. The 
good-deal pricing bounds enable us to avoid unreasonable prices and also to 
examine the price sensitivity to changes in the market price of risk. 

Sharpe Ratios are also used to price options 



In Bayraktar and Young 2008 



in incomplete markets. However, the perspective is that of an individual seller 
of one option, rather than that of the entire market. The seller of an option 
decides the option price via his own risk preferences, as expressed by his own 
chosen Sharpe Ratio. In other words, the seller of the option chooses the risk- 
neutral martingale measure under which he prices the option. It is shown in 
Bayraktar and Young 2008 that the upper and lower good-deal bounds of 
Cochrane and Saa Requejo 2000 can be obtained; in that case, the seller's 
chosen risk-neutral martingale measure coincides with the martingale measure 
which gives the upper good-deal bound. The lower good-deal bound is obtained 



Bayraktar and Young 2008 by considering the buyer of the option. 



A utility-based approach to the good-deal bound idea is found in Cerny 



2003 , and extended in Kloppel and Schweizer 2007 . An alternative approach 



based on the gain-loss ratio, which is the expectation of an asset's positive 
excess payoffs divided by the expectation of its negative excess payoffs, is found 



Bernardo and Ledoit 2000 



The aim of this paper is to apply the good-deal bound idea to the pricing of 
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derivatives in a regime-switching diffusion market. The paper is structured as 
follows. Section [2] details the regime-switching market model. In Section [3] we 
identify the set of equivalent martingale measures. In Section[4]thc Sharpe Ratio 
of an arbitrary asset in the market is defined and we state the extended Hansen- 
Jagannathan bound. The definitions of the upper and lower good-deal bounds 
on the price of a derivative are in Section [5j The stochastic control approach 
that we use to find them is detailed in Section [6j The minimal martingale 
measure, which we consider to be a benchmark pricing measure, is given in 
Section [Tj A numerical example illustrating the upper and lower good-deal 
bounds of a European call option of various maturities is in Section |8] We 
also examine the stability of the good-deal bounds when we change the market 
model's parameters. Finally, we conclude in Section [9] with some remarks. 

2 Market model 

We consider a regime-switching diffusion market model in which there is one 
risk-free asset traded asset and a finite number N of traded assets. An example 
of a risk-free asset is a bank account and typical examples of risky assets are 
equities, bonds or a pooled fund. The full mathematical description of the 
market is given below. 

2.1 Description of the market model 

We consider a continuous-time financial market model on a complete probability 
space (n, T, P) where all investment takes place over a finite time horizon [0, T], 
for a fixed T G (0,00). The probability space carries both a Markov chain a 
and an A'^-dimensional standard Brownian motion W = {Wi, . . . , Wn)^ , where 
we use to denote the transpose of A. 

The information available to the investors in the market at time t is the 
history of the Markov chain and Brownian motion up to and including time t. 
Mathematically, this is represented by the filtration 

Tt:^a{{a{s),W{s)),se[0,t]}vN{F), Vt e [0,T], (2.1) 

where A/'(P) denotes the collection of all P-nuU events in the probability space 
(fi, J",P). We assume that T = Tt- 

Remark 2.1. As a mathematical consequence of being defined on the same fil- 
tered probability space (51, J^, P, {J't}), the Markov chain and the Brownian 
motion are independent processes. Relating these processes to economic re- 
ality, we might think of the Brownian motion as modeling short-term, micro- 
economic changes in the market, whereas the Markov chain models long-term 
macro-economic changes. With this interpretation, the implicit assumption in 
our model that these economic changes are independent is a reasonable approx- 
imation to reality. For practical implementation, this means that the number 
and specification of the market regimes should be chosen to reflect this inter- 
pretation. 

Remark 2.2. In our model, an investor knows what regime the market is in at 
each time t. In reality, the market regime is unlikely to be known with certainty, 
although it could be estimated from market data. This is an important point 
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to note, since the prices of assets in the market are dependent on the initial 
market regime. 

The market is subject to regime-switching, as modelled by the continuous- 
time Markov chain a which takes values in a finite state space / = {!,..., D}, 
for some integer D > 2. For example, suppose we wish to model a market 
in which there are only two regimes: a bull market regime and a bear market 
regime. We set D — 2 and we might identify a{t) = 1 as corresponding to 
the market being in the bull market regime at time t. We would then identify 
a{t) = 2 as corresponding to the market being in the bear market regime at 
time t. 

We assume that the Markov chain starts in a fixed state io € /, so that 
a(0) = io, a.s. The Markov chain has a generator G, which is a. D x D matrix 
G = [9ij)?.j=i with the properties > 0, for aU i ^ j and gu = -J2j^i9ij- 
The interpretation of the off-diagonal clement gij of the generator matrix is as 
the instantaneous rate of transition from state i to state j. To avoid states 
where there are no transitions into or out of, we assume that ga < for each 
state i. 

Associated with each pair of distinct states (i, j) in the state space of the 
Markov chain is a point process, or counting process, 

5] l{a(._)=^}lMs)=,}, Vie[o,r], (2.2) 

0<s<t 

where 1 denotes the zero-one indicator function. The process Nij{t) counts the 
number of jumps that the Markov chain a has made from state i to state j up 
to time t. Define the intensity process 

^jii) '■= 9ij'^{a(t-)=i}- (2-3) 

If we compensate Nij(t) by Kjis) ds, then the resulting process 

M,,{t) -.^ N.,,{t) - f\,,{s)ds (2.4) 

^0 



is a martingale (see Rogers and Williams| [2006 Lemma IV. 21. 12]). We refer to 



the set of martingales [Mij] i, j G I ,i ^ j} as the V -martingales of a. They are 
mutually orthogonal, purely discontinuous, square-integrable martingales which 
are null at the origin. 

We consider a financial market that is built upon a finite number N of 
traded assets, which we call risky assets, and a risk-free asset. The risk-free rate 
of return in the market is denoted by the scalar stochastic process r and the 
risk-free asset's price process 5*0 — {5o(<),i e [0,T]} is governed by 

^^=r(Odi, Vie[0,T], 5o(0) = l. (2.5) 

The mean rate of return of the nth risky asset is denoted by the scalar stochastic 
process h„ and the volatility process of the nth risky asset is denoted by the 
A^-dimensional stochastic process <t„ = ((t„i, • • • , CTnjv)^. The price process 
Sn = {Sn{t),t E [0, T]} of the nth risky asset is then given by 

^^=5„(t)dt + a^(0dW(0, Vie[0,T], (2.6) 
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with the initial value S'„(0) being a fixed, strictly positive constant in K. 
Assumption 2.3. The market parameters r, h = (5i, . . . , ^Ar)^ and tr 



are sufBcicntly regular to allow for the existence of a unique 
strong solution to (2.51 and (2.6). Furthermore, the volatility process a is non- 
singular. 



3 Martingale measures 

3.1 Equivalent martingale measure 

From the fundamental theorem of arbitrage-free pricing, it is known that exis- 
tence of an equivalent martingale measure ("EMM") is equivalent to absence 
of arbitrage in the market. Furthermore, the market is complete (in the sense 
that all claims can be replicated) if and only if the EMM is unique. In our 
financial market model, while there is no arbitrage, the market is incomplete. 
This means that while EMMs exist, there is no unique one. This has imme- 
diate consequences for the valuation of contingent claims using our model, for 
example valuing European call options. We can price a European call option 
by the usual risk-neutral pricing formula. However, as there are infinitely many 
EMMs, we obtain a range of prices rather than a unique price. The good-deal 
bound approach is a means of narrowing the range of prices, which can be too 
wide to be useful in practice. The essential idea is to exclude those EMMs which 
imply a Sharpe Ratio that is too high. 



3.1.1 The Girsanov kernel process and the Girsanov Theorem 



Suppose we are given a probability measure Q on {D,,J-) which is equivalent 
to the (real-world) probability measure P. We define the likelihood process 
corresponding to the measure Q in the usual way as 



L{t) := E 



Tt], Vte[0,T]. 



We can assume that L is a positive {J'tj-martingale under the measure P (see 
Rogers and Wilhamsl [2006| Theorem IV.17.1]) with i(0) = 1, P-a.s. Recalling 
that the filtration is generated by both the Brownian motion W and the Markov 
chain a, we can apply an appropriate martingale representation theorem (for 
example, see Elliott 1976 Theorem 5.1]) to obtain predictable and suitably 
integrable stochastic processes (h,Tj), for h = {hi, ... , Hn)^ and rj :— {77^-; i,j — 
l,...,D,j ^ i}, satisfying 



dL{t) 



D D 



(t) dw(i) "^^1 d^'^^j- ' ^ ^] • 



(3.1) 



In order that the measure Q is non-negative, the process r\ must satisfy 

77,,(<)>-i, Vj^i, vte[o,r]. 

We call (h,J7) a Girsanov kernel process. As a consequence of the Girsanov 



theorem (for example, see Protter 2005 Theorem 40, page 135 
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we have 

dW{t) = h^{t) dt + dW^(t), (3.2) 
where, by Levy's Theorem, W'^ is a Q-Brownian motion; and 
the process 



Ml'^{t):^N,,{t)- il + rj,,{s))K,{s)ds, Vte[0,r], (3.3) 
Jo 

is a Q-martingale, for each j ^ i. We can interpret (1 + rjij{t)) Xij{t) as 
the intensity of the point process Nij{t) under the mcasm'e Q. 

The set of martingales {Mj^;i,j £ I,j ^ i\ arc the Q-martingales of a. They 
are mutuaUy orthogonal, purely discontinuous martingales which are null at the 
origin. Their integrability depends on the integrability of rjijit). Furthermore, 



substituting from (2.4 1 into ( |3.3[ ), we find 

dM,j(t) = v^J{t)>^^At) + dM%t), (3.4) 



which is analogous to (3.2) 



Remark 3.1. While a retains the Markov property under the measure Q (this 
can be shown using martingale problems, for example see Ethier and Kurtz] 



1986 Theorem 4.4.1]), it is not generally a Markov chain. This is because 
the intensity (1 + 77^ (t)) Xij{t) of the point process under the measure Q is not 
generally deterministic. 

Remark 3.2. The condition that riij{t) > —1 is to ensure that the measure 
Q is non- negative. However, if riij{t) = — 1 then P and Q are not necessarily 
equivalent which means that we can have arbitrage. As discussed in |Bjork and| 



Slinko 2006, Remark 3.3], to avoid any arbitrage possibility we can replace the 



inequality riij{t) > — 1 by r]ij{t) > — 1 + e, for some fixed < e <C 1 or we 
can regard any good-deal bounds derived with the constraint r]ij{t) > — 1 as 
open intervals of good-deal bounds. We choose the latter alternative since it is 
mathematically more convenient. 

3.1.2 The martingale condition and admissible Girsanov kernel pro- 
cesses 

Given a suitable process (h, 77) , we can generate a corresponding measure Q by 
using (3.1 1 to define the likelihood process L and then constructing the measure 
Qby 

^^m, on Ft. (3.5) 

Let Q be the measure generated by the Girsanov kernel process (h, jy). Consider 
an arbitrary asset in the market, with price process 11 = {n(i); t e [0, T]}. Note 
that this asset is not restricted to the traded risky assets or risk- free asset, 
but it could be any derivative or self-financing strategy based on them and the 
Markov chain a. If we price this asset using a risk-neutral measure Q, then the 
discounted price process is an {J^t}-martingale under the measure Q. As the 
filtration {Tt} is generated by both the Brownian motion and the Markov chain 



(recall (2.1)), then using a suitable martingale representation theorem (such as 
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Elliott 1976 Theorem 5.1]) this {J^t}-martingale can be expressed as the sum 
of a stochastic integral with respect to the Brownian motion and a stochastic 
integral with respect to the Q-martingales of the Markov chain. If we use the 
Girsanov theorem to obtain the P-dynamics of the price process, we still have 
a term involving the martingales of the Markov chain. This is the reason why 
the P-dynamics of the price process 11 are of the form 



dn(t) 
n(E) 



fi^it)dt+{a^{t)y dW{t) 



D D 

i=i 



(3.6) 



The processes fi^, ~ (o-p, . . . , cr]^)^ and are suitably integrable 

and measurable with the condition, in order to avoid negative asset prices, that 
^ —1 for each j ^ i. Note that if the asset is not the traded asset then the 
processes //^, and depend on the choice of the risk-neutral measure 

through the corresponding Girsanov kernel process. 



Apply (3.2) and (3.4 1 to (3.61 to obtain the price dynamics 11 of the arbi- 



trarily chosen asset under the measure 



dn(t) 
n(M 



D D 



f,^{t) + h^{t)a^{t) Kj{t) 
i=i j=i, 



dt 



(3.7) 



The measure Q is a martingale measure if and only if the local rate of return 
of the asset under the measure Q equals the risk-free rate of return r. Thus we 
obtain the following martingale condition. 

Proposition 3.3. Martingale condition The measure Q generated by the Gir- 
sanov kernel process (h,?j) is a martingale measure if and only if 

V^,{t)>^l, yj^i, (3.8) 
and for any asset in the market whose price process H has P-dynamics given by 



(3.6), we have 



D D 



r(t) = M"(t) + h^(t)an(t)+^^7n(^)^^^.(^)^^^.(^)^ Vie [0,r]. (3.9) 

We refer to a Girsanov kernel process (h, rf) for which the generated measure 
is a martingale measure as an admissible Girsanov kernel process. 



Remark 3.4. From (3.9) we have the following economic interpretation of an 



admissible Girsanov kernel process (hjjj): the process — h is the market price 
of diffusion risk and —rjij is the market price of jump risk, for a jump in the 
Markov chain from state i to state j. 
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Suppose we are given a Girsanov kernel process (h, 77) for which the generated 
measure Q is a martingale measure. The price dynamics under P of the nth 
underlying risky stock are as in (2.6), that is 

^Ml^b^(^t)dt + al{t)dW{t), Vie[0,T]. 

By Proposition |3.3[ we must have that 

r{t) ^bn{t)+h^{t)crn{t), Vte[0,r], for = 1, . . . , TV. 

This means that the market price of diffusion risk — h is determined by the price 
dynamics of the underlying risky assets, with the solution given by 

where 1 £ has all entries equal to one. However, as there is no traded asset 
in the market which is based on the Markov chain, we cannot say anything 
about the market price of jump risk —77^ . 



4 The Sharpe Ratio and a Hansen-Jagannathan 
Bound 

4.1 The Sharpe Ratio of an arbitrary asset 

We define a Sharpe Ratio process for an arbitrarily chosen asset, with P-dynamics 
as in (3.6 1. Broadly, the Sharpe Ratio is the excess return above the risk-free 
rate of the asset per unit of risk. We make this definition precise in our model. 
Define a volatility process v for the asset by 



d(n,n)(t) = n2(t_)i/2(i)di, 



(4.1) 



where (•, •) is the angle-bracket process. Substituting for 11 from (3.6) and using 
II- II to denote the usual Euclidean norm, we obtain 



/ 



d(n,n>(<)-n2(L) 



D D 



\ 



V 



=1 j=i, 



dt. 



(4.2) 



Comparing (4.1) and (4.2), we see that the squared volatility process satisfies 
.^(t) = ||an(0|P + f:f:i75(0pA„-(t). 

Recalling that the state space of the Markov chain a is denoted by / = {1, . . . , Z?} 
and the intensity process Xij{t) is given by (2.3), define the norm ||-||A(t) in the 
Hilbert space L'^{I x /,A(i)) by 

D D 



i=l 3 = 1, 
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Then we can write 

i.'it)^\\a^it)r + \\l''mlity 
Defining the Hilbert space 

H :=M^ X L2(/x /,A(i)), (4.3) 

and denoting by W-Wn the norm in the Hilbert space H, we can also express the 
volatility process as 

u{t)^\\{a^{t),j^{t))\\n. (4.4) 

Finally, we are in a position to define the Sharpe Ratio process {SR) for the 
arbitrarily-chosen asset. As is the local mean rate of return of the asset 
under the measure P, 

iSRW, (4^5) 

The Sharpe Ratio process depends on the chosen asset's price process. However, 
we seek a bound that applies to all assets' Sharpe Ratio processes. To do this, 
we use the extended Hansen- Jagannathan inequality, which is derived in [Bjork] 



and Slinko 2006 and is an extended version of the inequality introduced by 



Hansen and Jagannathan 1991 . 



4.2 An extended Hansen-Jagannathan Bound 

Lemma 4.1 (An extended Hansen-Jagannathan Bound). Recall the Hilbert 



space % in {4-3). For every admissible Girsanov kernel process (h,?;) and for 



any asset in the market whose price process H has ¥ -dynamics given by (3.6) 



and, consequently, whose Sharpe Ratio process {SR) is given by (4--5), the fol- 
lowing inequality holds. 

\{SR){t)\<\\{Ht)Mmn, 

that is 

D D 

|(^i?)(t)p<||h(t)|p + ^^|,y,^.(t)pA,,-(i). (4.6) 



Proof. The proof follows that of Bjork and Slinko 2006 Theorem A.l] and is 



therefore omitted. □ 
From Lemma |4.1[ we see that we can bound the Sharpe Ratios of all assets 



in the market by bounding the right-hand side of (4.6) by a constant. 



5 The general problem 



On the market detailed in Subsection |2.1[ we consider the valuation of a general 
contingent claim. To apply the good-deal bound idea, suppose we are given a 
contingent claim Z of the form 

Z:=a>(S(r),a(T)), (5.1) 
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for a deterministic, measurable function <I>, where S = {Si, . . . , Sn)'^ is the 
vector of the risky assets' price processes. As there is no unique martingale 
measure in the market, there is no unique price for the contingent claim. Rather 
than choosing one particular martingale measure to price the contingent claim, 
we seek instead to find a reasonable range of prices by excluding those martingale 
measures which imply Sharpe Ratios which are too high. 



5.1 The good-deal bound 

The key idea is that to restrict the set of martingale measures by way of the 
Sharpe Ratio, we use the Hansen-Jagannathan bound. Rather than bounding 



the Sharpe Ratios directly, we bound the right-hand side of (4.6 1 by a constant. 
We call the constant a good- deal bound. 

Condition 5.1. There exists Bq ^M. such that 

Bo = sup ||h(i)||2, a.s. 
te[o,T] 

Definition 5.2. A good-deal bound is a constant B > Bq. 

Remark 5.3. A chosen good-deal bound B bounds the Sharpe Ratio process 
(SR) of any asset in the market as follows: 

D D 

\{SR){t)\^<\\h{t)r + J2T.\'^^At)\'^vit)<B. (5.2) 

In other words, |(S'i?)(i)| < ^J~B. The economic interpretation is that, under the 
good-deal bound approach, \IB and —^/B are the highest and lowest achievable 
instantaneous Sharpe Ratio in the market, respectively. However, in the regime- 



switching diffusion market, we see from (5.2) that the good-deal bound B is 
really a bound on the price —rjij of regime change risk, since the price — h of 
diffusion risk is determined by the traded assets. 



5.2 The good-deal bound price processes 

We consider the problem of finding the upper and lower good-deal bounds on 
the range of possible prices of the contingent claim Z. 

Definition 5.4. Suppose we are given a good-deal bound B. The upper good- 
deal price process V for the bound B is the optimal value process for the control 
problem 

(5.3) 



supE^ ('e'-''*^'^(")^"$(S(T),a(r)) t\ 

(h,»)) V / 



(h,»)) 

where the predictable processes (h,Tj) are subject to the constraints 

Ht) = -{o''{t))-\h{t)-r{t)l), (5.4) 

nij{t)>~l, for i,j = 1, . . . ,D, j 7^ (5.5) 
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and 



D D 



(5.6) 



t=ij=i, 



for all t e [o,r]. 

Definition 5.5. The lower good-deal price process V is defined as in Definition 



5.4 except that "sup" in (5.3) is replaced by "inf" 



Remark 5.6. The risk-neutral valuation formula in (5.3 1 implies that the local 



rate of return of the price process corresponding to the contingent claim Z = 
^{S{T),a{T)) equals the risk- free rate r under the measure Q. The equality 
constraint ( |5.4[ ) ensures that h is consistent with the market price of jump risk. 
Together with the constraint (5.5), these ensure that the measure Q generated 

Note that 



by (h, T/) is a martingale measure, as in Proposition 3.3 



due to the 
generated by 



constant bound on (h,7j) in the constraint (5.6), the measure 
(h,»7) is a martingale measure, and not just a local martingale measure. 

Remark 5.7. The constraint (5.6) arises from the good-deal bound. It is ob- 



tained by bounding the right-hand side of (4.6). 

Remark 5.8. The only unknown in the constraints (5.4)-(5.6 ) is the market price 
of jump risk rjijit). If we obtain wide good-deal pricing bounds for a derivative 
then this tells us that the choice of the market price of jump risk riij[t) has a 
large impact on the derivative's price. Thus wide pricing bounds are a signal 
that we should explore additional ways of further restricting the possible values 



of the market price of jump risk rjij(t). This point is also made in Cochrane and 



Saa Requejo 2000 



The goal is to calculate the upper and lower good-deal bound price processes, 
which are what we consider to be reasonable bounds on the possible prices of the 
contingent claim Z. To calculate them, we use a stochastic control approach. 



6 Stochastic control approach 

To ensure that the Markovian structure is preserved under the martingale mea- 
sure Q, we need the following condition. 



Condition 6.1. The maximum in (5.3) is taken over Girsanov kernel processes 
(h, r}) of the form 

h(t) = h(t,S(t),a(i_)) and r^^j{t) = iq,j{t,^{t),a{t^)), Vj ^ z, 

and -quit) = 0, for all t e [0,T]. 



Remark 6.2. We note from the constraint (5.4) that the process h is completely 
determined by the market parameters r{t), b(i) and (j{t). This means that 
the requirement h(t) = h(t, S(t),Q;(t_)) is really a requirement that the market 
parameters are of the form 

r(t) = r(i,S(i),a(t_)), b(i) = b(t, S(t), a(i_)) and (T(t) = <T(t, S(t), a(i_)). 
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6.1 The good-deal functions 



Under Condition 6.1 the optimal expected value in (5.3) can be written as 



X / -> 



V(t,S{t), a{t-)) where the deterministic mapping V : [0,r] 
is known as the optimal value function. From general dynamic programming 
theory (for example, see Bjork 2009 Chapter 19]), the optimal value function 
satisfies the following Hamilton-Jacobi-Bellman equation 



dV 
~dt 



^up |a('''''V| -ry = 

y(T,x,z) = $(x,z), 



(6.1) 



where the supremum in (6.1) is subject to the constraints (5.4) - (5.6). An 
application of Ito's formula (for example, see Protter 2005 Theorem V.18, 



page 278]) shows that the infinitesimal operator A*^'^'''-' is given by 
Af*^-" V(t, X, z) 



^ dV 1 ^ 

(t,x, i) ^x„^^(t,x, i) + - ^ tT^(t,x, i)a„(t,x,i)a;„ 



n=l 



(i,x,i) 



(6.2) 



for all {t, X, i) G [0, T] 



X /. 



Definition 6.3. Given a good-deal bound -B, the upper good-deal Junction for 
the bound B is the solution to the following boundary value problem 



ay 



it,x,i) + sup \A^^'>^V{t.x,i)] -^r{t.x,i)V{t,x,i) = (6.3) 

(h,7j) J 

y(r,x,i) = <i>(x,j). 



where A^'^'''^ is given by (6.2) and the supremum is taken over all functions 
(h,r;) subject to Condition 6.1 and satisfying 



h{t,x,i) = - (o-^(t,x,i)) ^ ih{t,x,i) - r{t,x,i)l) , 
Vij{t,^) > -I7 for j = ■ ■ ■ ,D, j ^ i, 



and 



D 



\hit,X,^)\\^+J29^J\V^J{t,^)\''<B, 



(6.4) 
(6.5) 

(6.6) 



3 = 1, 



for all (t,x,i) e [0,T] x x /. We denote the solution to (6.3) by V^pp" 



Definition 6.4. The lower good- deal function is the solution to (6.3) but with 



the supr emu m replaced by an infimum, subject to Condition 6.1 
straints K^ - K^. We denote this solution by yi°^<=^ 



and the con- 
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Rather than attempting to solve the partial integro-differential equation 
("PIDE") of (6.3 1 directly, we reduce it to two deterministic problems which 
we solve for each fixed triple {t, x, i) S [0, T] x x /. Moreover, as h is 
completely determined by (6.4 1, we need to solve only for the optimal rj. 

Therefore, given h satisfying (6.4), we do the following. 

1. Solve the static optimization problem of finding the optimal fj in 



sup 

(h,»j) 



{A(h.'»)v^(t,x,z)} 



subject to the constraints (6.5) and (6.6). 
2. Using the optimal fj found above, solve the PIDE 



dV 
~dt 





$(x,i) 



(6.7) 



We consider in more detail how to solve the static optimization problem. To 
solve the PIDE, we can use numerical methods. A concrete example of this, 
where we find the good-deal bounds for a European call option, is given in 
Section [H 



6.2 The static optimization problem 

As we have seen above, the static optimization problem associated with the 
upper good-deal function of Definition 6.3 is to find for each triple (i,x, z) e 
[0, T] X X / the optimal i] that attains the supremum of 

N 



dV 1 ^ 



n=l 



n.m—l 



(t,x,i) 



(6.9) 



subject to the constraints 



(t,x)>-l, VjV* and ^5y|?7,,(t,x)p <B-||h(t,x,i)|l', (6.10) 



with h(t,x, i) given by (6.4). 

The static optimization problem associated with the lower good-deal func- 
tion of Definition 6.4 is as for the upper good-deal function but taking the 
infimum of (6.9), rather than the supremum. 
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As the only term in (6.9 ) which involves 77 is the last one, we can equivalently 



consider the problem of finding the optimal 77 which attains the supremum of 



D 



(6.11) 



subject to the constraints in (6.10). This is a linear optimization problem with 



both linear and quadratic constraints. We consider how the complexity of this 
problem increases as the number of states D of the Markov chain increases. 



6.2.1 Markov chain with tv^ro states 

When there are only two states of the Markov chain, the solution of the linear 
optimization problem (6.11) subject to the constraints (6.10) is very simple 
indeed and can be obtained by considering the sign of V{t,Ti,j) — V{t,x,i) in 
(|6lT). 



Lemma 6.5. For a 2-state Markov chain, fix [t, x, i) e [0, T] x M. x {1, 2} and 

define 



B(i,x,i) 



B - ||h(t,x,i) 

-gii 



1/2 



Then for each j ^ i, the solution to the static optimization problem associated 
with the upper good- deal function of Definition \6.^is 



B(^,x,^) 



l,S(t,x,i) 



ifV{t,^,j)-V{t,^,i)>Q 
tfVit,^,j)-V{t,^,i)<0, 



and the solution to the static optimization problem associated with the lower 
good- deal function of Definition\6.4\ is 



l,B(^,x,^) */y(t,x,j)-y(<,x,z) >o 



B{t,^,i) 



i/y(t,x,j)-y(t,x,z)<o. 



^■r(^,x) = 

Thus the solutions 7y'°™(t,x) and f^J^P^(t,x) depend on the value function 



V. This means that the numerical solution to the PIDE (6.7|-(6.8| involves 



checking at each node of the discretized state space the relative values of the 
value function in order to choose the appropriate solution. 



6.2.2 Markov chain with three or more states 

For a Markov chain with more than two states, the solution becomes more 
complicated because the number of constraints increases. If there are D states 
then for each fixed triple {t, x, i) S [0, T] x x / there are D — 1 variables 
{riij{t,'x) : j = l,...,D,j i} to find, each of which is subject to a lower 
and upper inequality constraint. Hence there are 2^~^ potential solutions, 
depending on which of the lower and upper constraints is binding. To see how 
the complexity increases, we consider a Markov chain with three states. For 
each X S K, denote by sgn(x) the sign of x. 
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Lemma 6.6. For a 3-state Markov chain, fix (t,x,i) G [0,r] x x {1,2,3} 
and define for each j e {1, 2, 3}, j ^ i, 



Then for each j =/= i and k ^ i, k ^ j. the solution (77j"^''(t, x), ^"^''(t, x)) to 
the static optimization problem associated with the upper good-deal function of 
Definition \6.3\ is one of the following pairs: 

(-1, -1), (6.12) 

(sgn(y(t,x,jVt"(i,x,z))B(t,x,z,j), -1), (6.13) 
(-1, sgn {V{t, X, k) - V{t, x, i)) Bit, X, ^, fc)), (6.14) 



(yit,^,j)-Vit,^,t)) 



(l/(i,x,fc)-l/(t,x,i)) 



g-||h(^,x,»)|p \ 
Xti9re{V{t,^,i}-V{t,^,i)fJ ' 

/■ \ 1/2 ■ 

B-\\h{t,^,i)r \ 

X?=i9uiV{t,^,£)-V{t,K,{}f) 



(6.15) 



and the solution {fjl°'^{t,x),fjl'^^{t,K)) to the static optimization problem asso- 
ciated with the lower good- deal function of Definition \6.4\ is one of the following 
pairs: 

(-1, -1), (6.16) 
(-sgn(y(^,x,J)-"^'(^,x,^))B(^,x,^,J), -1), (6.17) 
(-1, - sgn {Vit, X, k) - Vit, X, i)) B{t, x, z, k)), (6.18) 

{ - [Vit, X, ,) - V{t, X, z)) f , B-\\Ht,^,^W _\ 

V ' "\Y.t,9u{V{t,^,^)-V{t,^,^)f) 

^^ / S- ||h(i,X,z)|l2 \ \ 

VEfci5»n^(i>x,^) - y(<,x,i)) y y 

Proof. Apply the Kuhn- Tucker method. □ 

Remark 6.7. The solutions (77,"/P(i, x), 7?^"PP(t, x)) and {t,yi),ff^ {t,i^)) de- 
pend on the value function V , just as in the two-state Markov chain case. How- 
ever, the difficulty involved in solving the PIDE (6.7)-(6.8) numerically has 
increased since there are four potential solutions which must be checked at each 
node of the discretized state space. As we increase the number of states in the 
Markov chain, the number of potential solutions to the static optimization prob- 
lem increases and hence the complexity involved in solving the PIDE increases 
too. 
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7 Minimal martingale measure 



Here we leave aside the good-deal bounds and consider the minimal martingale 
measure, which we consider as a benchmark measure for pricing any derivative 
in the market. 

Definition 7.1. The minimal martingale measure is the measure Q'°"^ gener- 
ated by (h""",rj™'"), where (h™'",?;™'") is the Girsanov kernel process which 
minimizes 

D D 



subject to the constraint r(i,x, i) = 6„(t,x, i) 4- h^(t, x, i)£r„(t, x, i) for n — 
l,...,iV. 

It is immediate that the minimal martingale measure Q™'" is generated by 

\y^^\t):^-{o'(t)y\h{t)-r{t)\) and r;--(i) 0, VjV 

for all t e [0,r], where 1 e has all entries equal to one. As r/J"'"(i) > -1, 
we have that (h""",^™'") is an admissible Girsanov kernel process. 

Remark 7.2. Under the measure Q""", the process a is a Markov chain with the 
same generator G = {gij)^j^i as under the measure P. In particular, this means 
that the m easur e Q™'" preserves the martingale property of the process Mij (t) 
defined by (2.4), so that the P-martingales of a are also its Q™'"-martingales. 

Notice that (h™'",rj""") minimizes the right-hand side of (4.6) over the set 
of admissible Girsanov kernel processes. Moreover, by Dcfinition |5.2[ any good- 
deal bound B satisfies B > Bq. This means 

B>Bo^ sup ||h(i)||2= sup ||(h-'"(i),»7""'"(t))|p. 
te[o,T] te[o,T] 

Thus (h™'",?7™™) is a Girsanov kernel process which satisfies the good-deal 
bound constraint in (6.6). 

Denote the solution to the PIDE 

fIV 

+ A^h'"'"'"""" V - rl/ = (7.1) 
y(T,x,i) = $(x,z) (7.2) 

by Then as (h""",?;™'") is a Girsanov kernel process which satisfies the 

good-deal bound constraint in (6.6), it is clear from this and Definitions 6.3 and 



6.4 that the following relation holds: 

-flower ^ -^min ^ -^uppc 



8 Numerical example 

Having applied the good-deal bound idea in a regime-switching diffusion market, 
the next question is: are they useful? We examine this question by calculating 
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Table 1: Market parameters 



Regime i 


r{i) 


b{^) 


a{{) 


1 


0.06 


0.15 


0.12 


2 


0.06 


-0.22 


0.26 



the upper and lower good-deal pricing bounds for a 1-year European call option 
in a market where there are two regimes. We calculate them for various initial 
stock prices and for various choices of the good-deal bound. Finally, we examine 
how the pricing bounds change as we change the generator of the Markov chain 
which drives the regime-switching. 

8.1 Market model 

Suppose that we are in a financial market setting of Section |2] with only two 
market regimes and one risky asset, so that iV = 1, and time is measured in 
years. Assume the values of the market parameters given in Table [l] and take 
the generator of the Markov chain to be 

^ ^ (gii gi2\ ^ f-0-5 0.5 

U21 922) V 5 -5 



These figures are based on the estimated parameters found in Hardy 2001 for 
a 2-state regime-switching model fitted to data from the S&P 500, an index of 
500 U.S. stocks. From the generator G, we see that the average time spent in 
regime 1 is 2 years and the average time spent in regime 2 is about 2.5 months. 

8.2 Calculation and implementation 

We wish to calculate the upper and lower good-deal pricing bounds for a Eu- 
ropean call option with maturity T = 1 and strike price K = 100. To do this, 
we choose a good-deal bound B and fix the initial market regime a(0) = io 
and initial stock price S{0). Then we calculate the upper and lower good-deal 
functions of Definitions |6 .3| and 6.4 In Section[6] we saw that this involved first 



• solving the associated static optimization problem, and then 



• numerically solving the FIDE (6.7l-(6.8) using the solution to the static 
optimization problem. 

We have already solved the static optimization problem for a 2-state regime- 



switching model, with the solution given by Lemma 6.5 Thus it remains to 
numerically solve the FIDE 

dV , ^ , ^ dV , , In,, nd'^V , 

— (t,x,i)+r{i}x — (t,x,i) + -a {i)x -^{t,x,i) - r{i)V{t,x,i) 

-gu{l + m,{t,x)){Vit,x,j)-V{t,x,t))^0 (8.1) 
V(T, x,i) — max[a; — 0] , 

for j ^ i, using the optimal values fjij which solve the static optimization 
problem. Denoting the solution to the above FIDE by T/"pp°'' for the upper good- 
deal pricing bound and by for the lower good-deal pricing bound, the 
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good-deal price range for the derivative is (y'°^°''(0), y"PP°''(0)). The intervals 
are open intervals due to the discussion in Remark |3.2| 

We implement the numerical solution of the PIDE using a fully implicit, 
finite-difference method, based on a grid which has the values 



At = 0.01, AS* = 0.5, Sn 



0, 



200, 



where At is the grid step-size in the time direction (measured in years), AS is 
the grid step-size in the stock price direction and [S'min, •S'max] is the grid range 
in the stock price direction. The grid range in the time direction is [0,T]. We 
use the boundary conditions 

ViO,t) = and 1/(^„,ax,i) = S^,^ - Ke-'^'^^''\ Vt G [0,T]. 



Note that by Definition 5.2 and using the figures in Table [T] to calculate h{l) 
and h{2), the good-deal bound B must satisfy 

B > msix[h^{l),h^{2)] = max[(-0.750)^ (1.077)^] = 1.160. 
8.3 Results 

We begin by fixing the good-deal bound B = 1.2 and calculating the upper and 
lower pricing bounds for a range of initial stock prices. The results are shown 
in Figure [Tj with Figure 1(a) and 1(b) corresponding to the market starting in 
regime 1 and 2, respectively. The middle line in each of the plots corresponds to 
the minimal martingale measure price, which is the benchmark price. The plots 
show that, for the choice B = 1.2, the good-deal pricing bounds are reasonably 
narrow and therefore they are potentially of practical use. Next we examine 
exactly how the pricing bounds change as we vary the good-deal bound B. We 
fix the initial stock price 5(0) = 100 and calculate the good-deal pricing bounds 
for various choices of the good-deal bound B. These results are shown in Figure 
[2] with Figure |2(a) and 2(b) corresponding to the market starting in regime 
1 and 2, respectively. Again, the minimal martingale measure prices are the 
horizontal lines in the middle of each plot. The plots in Figure [2] show that as 
we increase the good-deal bound B, we increase the upper good-deal pricing 
bound. However, while the lower good-deal pricing bound decreases in Figure 



2(b) it is constant in Figure 2(a) The reason is that, in this particular market. 



the solution to the static optimization problem for the lower good-deal function 
is always ?7i2^(t, a;) = — 1 when starting in regime 1, regardless of the value 
of the good-deal bound B. Setting fji2{t,x) :— —1 in the PIDE (8.1), we see 
immediately that the last term on the left-hand side vanishes and hence the 
PIDE reduces to the classical Black-Scholes formula for a European call option 
in a non-regime-switching market with market parameters r{l), 6(1) and cr(l). 



8.4 Stability of the good-deal pricing bounds 

We base the market parameters in Table [T] on figures found in Hardy 2001 



However, the analysis in Hardy 2001 gives a large standard error in the estima- 

This leads us to wonder what happens if 
Do we have stability 



tion of the Markov chain parameters, 
we have mis-estimated the generator of the Markov chain, 
of the good-deal pricing bounds? To examine this issue, we consider again a 
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2 1 1 < < 1 < < 1 < 1 1 

90 92 94 96 98 100 102 104 106 108 110 
Initial stock price 

(a) Starting in regime 1. 




2 I 1 1 1 1 1 1 1 1 1 1 

90 92 94 96 98 100 102 104 106 108 110 
Initial stock price 

(b) Starting in regime 2. 

Figure 1: The upper and lower good-deal pricing bounds plotted against the 
initial stock price for a 1-year European call option with strike price K = 100 
and good-deal bound B = 1.2. The upper curve on each plot is the upper good- 
deal pricing bound, the lower curve is the lower good-deal pricing bound and 
the middle line is the minimal martingale measure price. The top plot assumes 
that the market is in regime 1 at time and the bottom plot assumes that the 
market is in regime 2 at time 0. 



19 



11.5 



11 - 



10.5 



0- 9.5 - 



8.5 



1.6 



1.8 2 2.2 2.4 

Good-deal bound 



2.6 



2.8 



(a) Starting in regime 1. 



11.5r 



11 - 



10.5 



°- 9.5 ■ 



i.5 - 



1.2 



1.4 



1.6 



1.8 



2 2.2 
Good-deal bound 



2.4 



2.6 



2.8 



(b) Starting in regime 2. 



Figure 2: The upper and lower good-deal pricing bounds plotted against the 
good-deal bound for a 1-year European call option with strike price K = 100. 
The initial stock price is 5(0) = 100. The top plot assumes that the market 
is in regime 1 at time and the bottom plot assumes that the market is in 
regime 2 at time 0. On both plots, the minimal martingale measure price is the 
horizontal line in the middle. 
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Table 2: Diagonal elements of the generator G of the Markov chain. 





Model 1 


Model 2 


Model 3 


-511 


0.5 


0.333 


0.667 


-522 


5 


6 


6 


Avg. time in regime 1 


2.000 


3.000 


1.500 


Avg. time in regime 2 


0.200 


0.167 


0.167 



1-year European call option with strike price K = 100. We find the good-deal 
pricing bounds for this option for three different models, assuming that the price 
of the risky stock at time is S{0) = 100. In each of the models, the market pa- 
rameters r, b and a are as in Table [T] but the diagonal elements of the generator 
G of the Markov chain are given in Table [2j Note that Model 1 corresponds to 
the model described in Subsection |8.1[ The results are shown in Figure |3j with 
Figure 3(a) and |3(b)] corresponding to the market starting in regime 1 and 2, 
respectively. The middle horizontal lines in the plots correspond to the minimal 
martingale measure prices. 

The plots show that the good-deal pricing bounds are sensitive to the choice 
of the generator G of the Markov chain. Roughly, the upper good-deal pricing 
bounds move in tandem with the minimal martingale measure prices. In Figure 
the lower good-deal pricing bounds behave similarly. However, in Figure 
the lower good-deal pricing bounds are nearly all constant. The explana- 



3(b) 



3(a) 



tion for the constant lower pricing bounds is as before: the solution to the static 
optimization problem for the lower good-deal function is Tf°2'{t,x) = — 1 when 
starting in regime 1 in these cases. In Figure [3 (a) | we see that lower good-deal 
pricing bound is slightly higher at the good-deal bound B = 1.2 because here 
the solution Tfi^{t, x) to the static optimization problem for the lower good-deal 
function is just above —1. 



9 Conclusion 



We have applied the good-deal bound idea of Cochrane and Saa Requejo 2000 



to a regime-switching market using the approach of Bjork and Slinko 2006 and 



illustrated it with a numerical example. The good-deal bound idea is a way to 
measure the uncertainty in the choice of the risk-neutral measure used to price 
derivatives. However, as our numerical example demonstrates, the good-deal 
pricing bounds change when the model changes, even though the good-deal 
bound remains constant. Thus the good-deal pricing bounds are sensitive to 
the choice of model. It would be interesting to do a wider investigation of 
the variation of the good-deal pricing bounds over a wider class of models for 
different derivatives. 

We have looked at pricing, but what is a "good-deal hedging" strategy? As 
Bjork and Slinko 2006 say in their conclusion, this is a highly challenging open 
problem. 
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Model 1 




8.5 - 



8 - 



1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 

Good-deal bound 



(a) Starting in regime 1. 



11.5 




8.5' ' ' ' ' ' ' ' ' ' 

1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 

Good-deal bound 

(b) Starting in regime 2. 

Figure 3: The upper and lower good-deal pricing bounds plotted against the 
good-deal bound for a 1-year European call option with strike price K := 100. 
The bounds are calculated for three market models with identical market pa- 
rameters (given in Table [ij but with different generators G of the Markov chain 
(with diagonal elements given in Table [2]). The top plot assumes that the mar- 
ket is in regime 1 at time and the bottom plot assumes that the market is in 
regime 2 at time 0. For each model, the minimal martingale measure price is 
the horizontal middle line. 
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